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基于 无 人 机 遥感 表 型 监测 的 芝麻 优质 种 质 资源 筛选 方法 


taun, E SE A BR, BAA, PAR, ERR, 
陈 建 福 , A H, Be 
(湖南 农业 大 学 农学 院 ,湖南 长 沙 410128) 
摘 m. 苍 麻 是 重要 的 纤维 作物 之 一 ， 由 于 土地 资源 紧缺 及 优良 品种 的 推广 应 用 等 原因 ， 茸 麻 遗 传 变异 和 


遗传 多 样 性 减少 ， 对 芝麻 种 质 资 源 多 样 性 调查 和 保护 的 需求 日 趋 加 大 。 基 于 无 人 机 遥感 的 作物 表 型 测量 方 
法 可 以 对 不 同 基 因 型 作物 的 生长 特性 进行 频繁 、 快 速 、 无 损 、 精 准 的 监测 ， 实 现 作 物种 质 资源 调查 ， 筛 选 
特异 优质 品种 。 o o CM 辅助 筛选 优势 芝麻 品种 ， 本 研究 提出 了 一 种 
基于 无 人 机 遥感 影像 的 芝麻 种 质 资源 表 型 监测 及 筛选 方法 。 首 先 ， 基 于 无 人 机 遥感 影像 ， 利 用 Pix4dmapper 
软件 生成 试验 区 的 数字 地 表 模 型 (Digital Surface Model, DSM) 和 正 射 影像 ; 然后， 对 蔚 麻 种 质 资源 关键 表 
型 参数 〈 株 高 、 株 数 、 叶 面积 指数 、 叶 片 叶绿素 含量 、 含 水 量 ) 进行 估 测 。 基 于 DSM 采 用 “差分 法 ”提取 
芝麻 株 高 ， 基 于 正 射 图 像 采 用 目标 检测 算法 提取 芝麻 株数 ， 采 用 机 器 学 习 方 法 佑 测 芝麻 叶 面 积 指 数 (Leaf 
Area Index, LAI), 、 叶 片 叶 绿 素 含 量 (SPAD 值 ) 、 含 水 量 ; 最 后 ， 根 据 提取 的 各 项 遥感 表 型 参数 ， 采 用 变异 
性 分 析 和 主 成 分 分 析 方法 对 芝麻 和 中 质 资源 进行 遗传 多 样 性 分 析 。 结果 表明 ，(1) 基于 无 人 机 遥感 的 芝麻 表 
型 估 测 效果 较 好 ， 株 高 的 拟 合 精度 为 0.93， 均 方 根 误差 为 5.65 em; SPAD 值 、 含 水 量 、LAI 的 拟 合 指标 分 别 
达到 0.66、0.79、0.74，RMSE 分 别 为 2.03、2.21、0.63; (2) 芝麻 种 质 资源 的 遥感 表 型 存在 较 大 差异 ，LAI、 
株 高 和 株数 的 估 测 值 变异 系数 分 别 达到 20.82%、24.61% 和 35.48%; (3) 利用 主 成 分 分 析 法 将 蔚 麻 种 质 资源 
的 遥感 表 型 聚 类 为 因子 1 CERES. LAT) 和 因子 2 (LAI、SPAD 值 )， 因 子 1 可 用 于 区 麻 种 质 资 源 结 构 特 征 评 
价 ， 因 子 2 可 以 作为 高 光 效 世 麻 资源 的 筛选 指标 。 本 研究 将 为 作物 种 质 资源 表 型 监测 和 育种 相关 分 析 提 供 
参考 。 
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土壤 修复 等 方面 得 到 广泛 应 用 。 中 国 是 芝麻 品种 
变异 类 型 和 芝麻 属 野 生 种 较 多 的 国家 ， 但 随 着 十 
地 资源 紧缺 以 及 优良 品种 的 推广 应 用 ， 这 一 宝贵 
的 作物 资源 出 现 严重 的 品种 单一 化 趋势 ， 导 致 艺 
麻 种 质 资源 面临 威胁 ， 对 芝麻 种 质 资源 多 样 性 调 
查 和 保护 的 需求 日 趋 加 大 。 

作物 种 质 资 源 中 列 含 着 丰富 的 遗传 多 样 性 ， 
并 且 在 长 期 环境 适应 及 人 工 选择 的 互 作 影响 下 ， 
形成 了 表 型 性 状 的 多 样 性 ， 如 农艺 性 状 、 产 量 性 
状 、 抗 性 性 状 、 品 质 性 状 等 ， 构 成 选 育 作 物 新 品 
种 的 物质 基础 ”“。 面 对 成 千 上 万 个 种 质 资 源 ， 表 
型 数据 的 测量 是 一 件 巨 大 繁 珊 的 工程 。 传 统 的 人 
工 测量 方法 不 仅 耗 时 耗 力 ， 并 且 受 环境 和 观察 
主观 因素 影响 较 大 ， 难 以 实现 多 地 点 、 多 时 序 的 
种 质 资源 表 型 测量 工作 。 高 通 量 表 型 技术 的 发 展 
为 作物 表 型 性 状 评估 提供 新 途径 7, d Marefat- 
zadeh 等 ”利用 半自动 室内 高 通 量 技术 对 来 自 不 
同 地 区 的 589 份 番茄 种 质 资源 进行 了 特征 分 析 ， 
证 明了 高 通 量 表 型 技术 在 作物 种 质 资源 遗传 多 样 
性 研究 中 的 应 用 潜力 。 

随 着 室内 高 通 量 表 型 技术 的 日 益 成 熟 ， 该 项 
技术 的 应 用 逐渐 由 受 控 环境 转向 田间 复杂 环境 ， 
越 来 越 多 的 研究 者 开始 关注 高 通 量 表 型 技术 在 田 
间作 物 生长 监测 上 的 应 用 。 田 间 表 型 技术 依托 无 
人 机 遥感 技术 、 光 谱 技 术 、 人 工 智 能 等 技术 的 发 
展 ， 在 作物 株 高 测量 ”、 品 种 识别 ARB 
Dor  、 产 量 估 测 ”“、 病 虫害 监测 ”等 方面 取得 
突破 。 如 Jiang 等 "" 基于 无 人 机 多 光谱 遥感 获取 
的 植被 指数 ， 采 用 随机 森林 回归 方法 估 测 蒙 麦 叶 
面积 指数 (Leaf Area Index, LAI) 和 叶绿素 含 
量 ， 佑 测 精 度 分 别 为 0.977~0.980 和 0.983— 
0.986。Weiss FI Baret "基于 无 人 机 遥感 系统 获 
取 葡 萄 树 树 冠 影像 ， 采 用 三 维 重 建 技术 提取 了 和 葡 
萄 树 高 度数 据 ， 其 精度 可 达到 88%。 Ganeva 
等 "基于 无 人 机 遥感 图 像 ， 在 育种 试验 中 精确 
提取 了 冬 硬 粒 小 麦 LAI、 吸 收 光 合 有 效 辐 射 分 
数 、 植 被 覆盖 分 数 、 叶 片 叶绿素 含量 (SPAD 
值 ) 、 冠 层 叶绿素 含量 、 分 蓝 和 产量 ， 并 分 析 了 


所 研究 基因 型 的 遗传 多 样 性 、 接 近 性 和 相似 性 。 
Tanger 等 ”人 研究 了 高 通 量 表 型 技术 用 于 识别 水 
稳重 要 性 状 遗 传 变 异 的 潜力 ， 发 现 高 通 量 表 型 技 
术 获 得 性 状 识别 的 数量 性 状 位 点 结果 与 人 工 观测 
PEAR AY QTL (Quantitative Trait Locus) 结果 相 
对 应 。 

当前 ， 基 于 无 人 机 遥感 的 作物 种 质 资源 田间 
表 型 高 通 量 监测 及 辅助 育种 研究 仍 处 于 初步 探索 
阶段 ， 以 往 研 究 往 往 受 限于 作物 品种 ， 或 侧重 对 
单一 表 型 信息 进行 监测 ， 缺 乏 对 作物 关键 表 型 的 
综合 考虑 。 因 此 ， 本 人 研究 以 更 多 类 型 的 芝麻 种 质 
资源 为 对 象 ， 根 据 芝 麻 种 质 资源 评价 中 常用 的 表 
型 考察 指标 ， 探 究 了 无 人 机 迁 感 技术 在 芝麻 种 质 
资源 多 样 性 调查 和 评价 中 的 可 行 性 ， 并 探索 利用 
遥感 技术 得 到 的 作物 表 型 分 析 结 果 筛 选 优质 品 
种 、 辅 助 作物 育种 的 可 行 性 。 主 要 工作 为 利用 无 
人 机 遥感 系统 获取 田间 芝麻 冠 层 影像 ， 在 此 基础 
上 采用 图 像 处 理 技术 、 机 器 学 习 技 术 提 取 臣 麻 种 
质 资源 主要 表 型 ( 株 高 、 株 数 、LAI、SPAD 值 、 
含水 量 ) ， 最 后 基于 遥感 表 型 对 芝麻 种 质 资源 进 
行 遗传 多 样 性 和 差异 研究 。 


2 材料 与 方法 


2.1 研究 区 域 


研究 区 域 位 于 湖南 省 长 沙市 美容 区 湖南 农业 
KE te bd PPE HE (28°11'02"N, 113°04'10"E), 
DX ASE RAT 154 0) Ras f (图 1)， 于 2017 年 
12 月 育苗 移 栽 ，2018 年 6 月 首次 收获 。 区 域内 共 
包含 1$4 个 小 区 ， 各 小 区 面积 为 3.6m， 小 区 内 
种 植 8 划 这 及 (2 行 X4 昔 )， 基 因 型 间距 离 为 
0.6 m， 排 水 沟 宽 为 0.5 m。 试 验 区 土壤 成 分 均一 
且 土 壤 肥 沃 ， 灌 溉 排水 便捷 ， 田 间 水 肥 管 理 一 
致 。 为 完成 后 期 影像 的 重建 ， 提 高 经 纬 、 高 程 位 
置 校准 的 精度 ， 试 验 地 布置 了 6 个 地 面 控制 点 。 


2.2 数据 采集 


于 2022 年 芝麻 苗 期 (3 月 15 日 、3 月 23 日 )、 
封 行 期 (3 月 29 日 、4 月 7 日 )、 旺 长 期 (4 月 12 
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A» X 
图 1 芝麻 种 质 资源 试验 小 区 位 置 及 分 布 图 
Fig.1 Location and distribution map of ramie germplasm 


resources test plot 


日 、4 月 20 日 ) 进行 了 6 次 数据 采集 ， 包 括 地 面 
测量 数据 和 无 人 机 影像 数据 。 为 保证 数据 的 时 效 
性 ， 地 面 数据 测量 与 无 人 机 航拍 在 同一 天 同步 
进行 。 
2.2.1 地面 数据 采集 

采集 了 作物 株 高 、SPAD 值 、LAI、 含 水 量 、 
株数 作为 地 面 数据 。 其 中 ， 株 高 利用 直 尺 测量 芝 
麻 植 株 底部 至 冠 层 项 端的 距离 ; SPAD 值 通过 日 
AX KONICA MINOLTA 公司 生产 的 SPAD-502 型 
叶绿素 仪 测量 ; LAI 采 用 LAI-2200 冠 层 分 析 仪 测 
量 。 叶 片 含 水 量 为 叶片 鲜 重 与 供 干 叶片 重量 的 差 
值 。 以 上 数据 均 通 过 随机 选取 各 小 区 10 株 芝麻 
求 取 平均 值 获取 。 于 2022 年 4 月 20 日 ， 通 过 人 
工 目 视 方法 获取 各 小 区 芝麻 植株 数量 。 
2.2.2 ”遥感 影像 采集 

采用 大 疆 悟 2 搭载 蝉 思 XSss 高 清 数码 相机 获 
取 芝 麻 冠 层 RGB 影像 ， 采 用 大 性 精灵 4 获取 芝麻 
冠 层 多 光谱 影像 。 

为 保证 太阳 辐射 稳定 和 充足 ， 飞 行 作业 均 选 
择 在 当地 时 间 12:00—14:00 进行 。 采 用 DJI-GS 


Pro 在 指定 区 域内 自动 生成 航线 ， 设 置 主 航向 重 
BE 80%, EZR ARB 70%, KTR 
为 20m， 云 台 俯 仰角 为 -90"， 相 机 曝光 模式 为 
自动 。 


2.3 遥感 影像 预 处 理 


2.3.1 遥感 影像 拼接 

采用 Pix4dmapper 软件 生成 研究 区 域 的 正 射 
影像 以 及 数字 地 表 模 型 (Digital Surface Model, 
DSM) 影像 (图 2)。 拼 接 过 程 中 ,为 获取 准确 
的 地 理 参考 ， 对 导入 地 面 控制 点 的 三 维 空间 位 置 
信息 进行 地 理 校正 。 


(a) 正 射影 像 (b) 数 字 地 表 模 型 


图 2 芝麻 种 质 资源 试验 区 域 正 射影 像 及 数字 地 表 模 型 


Fig. 2 Orthophoto image and digital surface model of ramie 


germplasm resources test area 


2.3.2 ”遥感 特征 值 提取 

从 无 人 机 遥感 影像 中 提取 的 遥感 特征 值 包括 
纹理 特征 值 、 光 谱 特 征 值 以 及 H,w (高 程 数 据 )。 
选用 灰 度 共生 和 矩阵 作为 纹理 特征 值 ， 植 被 指数 作 
为 光谱 特征 值 ， 植 被 指数 通过 对 不 同 波段 的 反射 
率 进 行 线性 或 非 线性 方式 组 合 而 成 “5 ;Hosw 通 
过 对 不 同时 期 获取 的 DSM 进行 差 值 运算 提取 。 
具体 计算 公式 见 表 1。 

通过 绘制 矩形 感 兴趣 区 域 (Area of Interest, 
AOD 提取 单位 小 区 的 平均 遥感 特征 值 。 采 用 
Arc GIS10.2 软件 中 的 栅 格 工具 ， 治 小 区 田 块 边 
缘 绘 制 AOI。 绘 制 时 ， 小 区 四 周边 缘 余 留 10% 以 
消除 边缘 效应 干扰 。 
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TT 
表 1 遥感 特征 值 的 计算 公式 
Table 1 The calculation formula of remote sensing index 
类 型 遥感 指数 计算 公式 
灰 度 共 生 和 矩阵 均值 (Mean) Mean = us PiN? (1) 
XJ EE RE ( Contrast) Contrast = > ， yl i-j) ? pi, j) (2) 
255 (Dissimilarity ) Dissimilarity = p» us "ur li - jl (3) 
[HIE (Homogeneity ) Homogeneity = 2. A : y p(ij) (4) 
纹理 特征 值 
能 量 (Energy) Energy = b» (5j) (5) 
相关 性 (Correlation ) Correlation = 220 n) U-» (6) 
fti — BY FE (Angular Second Moment, ASM) ASM = Spi jy (7) 
y= ‘ei Entropy ) Entropy = >), >) p (inj) log p (ij) (8) 
a 归 一 化 植被 指数 (Normalized Difference Vegetation Index, NDVI) NDVI = (NIR-R)/( NIR+R) (9) 
eo 绿色 归 一 化 差异 植被 指数 (Green Normalized Differential Vegetation 
GNDVI = (G-R)/(G+R) (10) 
= Index, GNDVI) 
= 土壤 调节 植被 指数 (Soil Adjusted Vegetation Index, SAVI) SAVI = 1.5(NIR-R)/(NIR+R+0.5) (11) 
" É 曾 强 植被 指数 (Enhanced Vegetation Index, EVI) EVI-2.5(NIR-R)/( NIR+6R-7.5B+1) (12) 
N 光谱 特征 值 g 
e 超 绿 指数 (Exceed Green Index , ExG) ExG 2 2X g-r-b (13) 
eo 超 红 指数 (Exceed Red Index, ExR ) ExR = l.4r-g (14) 
过 绿 红 指 数 (Exceed Green And Red Index, ExGR) ExGR = ExG-1.4R-G (15 
É 寸 绿 红 指数 ) 
N 可 见 大 气 抗 性 指数 (Visible Atmospheric Resistance Index, VARI) VARI = (g-r)/Cg*r-b) (16) 
> Hypon 提取 的 高 程 数据 (HDSM) Hpsu= DSMi-DTM (17) 
x 注 : 1 表示 像素 点 行 数 ; /为 像素 点 列 数 ; p (i, D) 表示 图 像 中 (i, j) BEKE; ox, oy Ap, (0) fep, O) 的 方差 ，N 为 行 或 列 
© 3k; NIR、R、B、G 分 别 表示 近 红 外 波段 、 红 波段 、 蓝 波段 、 绿 波段 的 平均 反射 率 ; 7 表示 红色 标准 值 ; g 表 示 绿 色 标 准 值 ; bp 表示 蓝 
c BARE, DSMi cc SY iM AS REZ DSM, DIMAS AA KA} HY DSM 
= 2.4 Pa ARE I ze 8 M Ml) 


2.4.1. ERRERA A 

在 前 期 研究 中 ,利用 2019 年 和 2020 年 获取 
的 芝麻 种 质 资源 冠 层 影像 ， 采 用 全 卷 积 一 阶段 目 
标 检 测算 法 (Fully Convolutional One-Stage Ob- 
ject Detection, FCOS) 构建 了 芝麻 植株 计数 模 
型 。FCOS 模型 采用 了 Resnet101 作为 backbone， 特 
征 金字 塔 (Feature Pyramid Networks, FPN) 用 


于 提取 影像 中 的 芝麻 植株 特征 ， 训 练 次 数 为 20 图 3 基于 全 卷 积 一 阶段 目标 检测 算法 的 

次 。 图 3 为 芝麻 植株 检测 结果 ， 蓝 框 标注 实际 植 芝麻 植株 检测 结果 

株 样 本 ， 红 框 标注 识别 样本 。 Fig. 3 Ramie plant counting results based on fully convolu- 
24.2 ”区 麻生 理 参数 监测 tional one-stage object detection(FCOS ) 


混合 各 时 期 采集 的 样本 数据 ， 按 7:3 比例 划 
分 为 训练 集 和 验证 集 。 以 往 研 究 表明 ， 作 物 生 理 
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生化 指标 含量 对 不 同 遥 感 特征 的 响应 存在 差异 ， 
如 作物 LAI 被 证 明 与 相对 活力 指数 (Relative 
Vigor Index，RVI) 、 转 换 型 土壤 调节 植被 指数 
(Transform Soil Adjusted Vegetation Index, 
TSAVI) 呈 线 性 相关 ,与 绿色 植被 指数 
(Green Vegetation Index，GVI) 、 垂 直 植 被 指数 
(Perpendicular Vegetation Index，PVI) 更 适合 
二 次 方程 拟 合 。 含 水 量 与 ExG 、ExR 等 植被 指数 
具有 极 显著 相关 性 ""。 本 人 研究 以 田间 实测 表 型 
数据 和 遥感 特征 值 为 输入 参数 ， 递 归 特 征 消 除 方 
法 筛选 敏感 特征 ， 采 用 4 种 广泛 使 用 的 机 器 学 习 
方法 : 线性 回归 (Linear Regression, LR), BAHL 
森林 回归 (Random Forest Regression, RF), x 
持 向 量 机 回归 (Support Vector Machines, 
SVM) 、 偏 最 小 二 乘 回归 (Partial Least Squares 
Regression Analysis, PLSR) 构建 芝麻 生理 参数 
估 测 模型 。 使 用 决定 系数 (Coefficient of Deter- 
mination, R’), JIJINI (Root Mean Square 
Error, RMSE) 作为 模型 评价 指标 ，R 越 高 ， 
RMSE 越 低 ， 表 明 模 型 拟 合 精度 越 高 。 


2.5 芝麻 种 质 资源 遗传 多 样 性 分 析 


传统 芝麻 种 质 资源 评价 方法 存在 费时 费力 、 
破坏 性 强 等 缺陷 ， 为 此 ， 本 研究 采用 遥感 监测 得 
到 的 表 型 数据 对 芝麻 种 质 资源 表 型 遗传 多 样 性 和 
遗传 变异 进行 分 析 ， 以 实现 作物 种 质 资源 的 高 效 
评价 与 第 选 。 变 异 系数 和 和 多样 性 系数 为 遗传 
多 样 性 提供 了 信息 支撑 ， 变 异 系数 越 高 ， 遗 传 多 
样 性 系数 越 高， 遗传 多 样 性 越 显著 。 主 成 分 分 析 
不 仅 能 够 提供 表 型 性 状 之 间 的 相关 性 信息 ， 还 能 


分 比 ，%; Sd 表示 标准 偏差 ， Mean 表示 平均 值 。 
3 结果 与 分 析 


3. 


-— 


eS DR PR Ja RE He EK eT 


基于 不 同 生 育 期 获取 的 地 面 实测 数据 对 芝麻 
种 质 资源 表 型 性 状 进行 统计 描述 ， 结 果 见 表 2。 
从 苗 期 到 旺 长 期 ， 苇 麻 种 质 资源 的 SPAD 值 和 含 
水 量 总 体 呈 现下 降 趋势 ， 这 一 变化 趋势 与 营 麻 生 
长 规律 吻合 。 叶 片 含水 量 是 反映 水 分 供应 与 蒸腾 
之 间 平 衡 关 系 的 重要 指标 ， 从 苗 期 到 旺 长 期 ， 气 
温和 逐渐 升 高 ， 干 旱 胁 迫 下 芝麻 叶片 凌 蓄 度 增 加 ， 
叶片 相对 含水 量 下 降 。 伴 随 作物 生长 ， 作 物 叶片 
数量 和 面积 的 改变 ， 使 得 LAI 处 于 连续 动态 变化 
之 中 。 本 人 研究 中 ， 从 苗 期 到 旺 长 期 ， 地 面 实际 测 
量 的 LAI 逐 渐 增 大 ， 这 是 由 于 随 着 麻 葵 生长 ， 出 
叶 总 数 逐 渐 增 加 ， 叶 片面 积 逐 渐变 大 ， 蔚 麻 行 间 
荫蔽 加 重 ， 因 此 LAI 也 随 之 持续 增长 AER PRI 
整体 表现 出 随 生 育 期 后 移 ， 逐 渐 增 高 的 趋势 。 

从 不 同时 期 芝麻 种 质 资源 存在 的 表 型 差异 来 
看 ,前 3 个 生育 期 内 ， 营 麻 叶 片 SPAD 值 总 体 在 
24.70~46.51 的 范围 ， 苗 期 SPAD 值 表现 出 明 品 
的 品种 差异 ， 变 异 范 围 为 23.50~46.51， 标 准 差 
达到 4.02， 这 表明 芝麻 种 质 资源 在 苗 期 光 效 利用 
上 存在 较 大 变异 ， 苗 期 可 能 是 影响 芝麻 有 机 质 最 
终 累 积 量 的 关键 时 期 。 随 着 芝麻 生长 ，SPAD 值 
在 品种 间 的 差异 减 小 ， 封 行 期 芝麻 种 质 资源 
SAPD 值 的 标准 差 为 3.39， 旺 长 期 为 3.36; 芝麻 
种 质 资源 含水 量 在 前 3 个 时 期 的 总 体 范围 为 
60.81% 一 95.84% ， 按 不 同时 期 芝麻 种 质 资 源 含 水 


够 根据 主要 成 分 的 个 体 性 状 值 提供 基因 型 定位 


量 标准 差 排序 ， 表 现 为 旺 长 期 > 封 行 期 > 苗 期 。 旺 


言 息 。 
以 上 分 析 均 在 SPSS 软件 中 实现 ， 变 异 系数 
和 多 样 性 系数 的 计算 公式 如 下 : 
CV=Sd/Mean X 100% (18) 
H'=> PiX InPi (19) 
其 中 ，CV 为 变异 系数 ; HH EEE RŽ 
Pi 为 某 表 型 性 状 第 i 级 内 个 体 数 占 总 个 体 数 的 百 


长 期 的 含水 量 差异 范围 为 60.81%~85.06%， 封 
行 期 芝麻 种 质 资源 的 含水 量 在 71.31% 一 93.49% 
的 范围 ， 苗 期 芝麻 种 质 资源 含水 量 范 围 为 
79.69%~-95.84%。 不 同 生 育 期 内 ， 芝 麻 种 质 资源 
表现 出 的 高 含水 量 差异 有 助 于 筛选 干旱 条 件 下 耐 
旱 性 强 的 品种 ; 前 3 个 生育 期 内 ， 芝 麻 种 质 资源 
LAI ŠRF 0.91~7.64 的 范围 ， 在 各 时 期 间 的 变 
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表 2 区 麻 种 质 资源 表 型 性 状 统计 描述 
Table 2 Statistical description of ramie phenotypic characters 
bp 生育 期 最 小 值 最 大 值 平均 值 标准 差 
苗 期 25.50 46.51 37.01 4.02 
. 封 行 期 24.70 44.52 36.76 3.39 
SPAD 值 
旺 长 期 26.80 44.60 35.49 3.36 
全 生育 期 24.70 46.51 36.55 3.66 
苗 期 79.69 95.84 87.73 2.36 
封 行 期 71.31 93.49 79.17 2.53 
含水 量 /% 
旺 长 期 60.81 85.06 76.96 3.07 
全 生育 期 60.81 95.86 80.23 4.89 
苗 期 0.91 5.03 2.77 0.76 
SH 1.05 7.25 3.60 1.03 
S ad 旺 长 期 2.24 7.64 4.64 0.91 
er 全 生育 期 0.91 7.64 3.65 1.19 
CO 苗 期 14.50 71.60 31.72 10.28 
T 封 行 期 21.70 178.90 58.07 20.38 
= TER 旺 长 49.60 252.10 92.49 22.42 
N | 全 生育 期 14.50 252.10 58.27 3051. 
e 株数 / 株 一 一 11 140 67 24 
co 
N 化 幅度 不 大 。 随 着 生育 期 的 后 移 ， 芝 麻 种 质 资源 RI 芝麻 种 质 资源 株 高 拟 合 模型 
E 间 株 高 差异 逐渐 变 大 ， 标 准 差 由 苗 期 的 10.28 变 Table 3 Plant height fitting model oframie 
S 化 到 旺 长 期 的 22.42。 germplasm resources 
d m 
~ A EER EE E EAM TEmo RAMA Sam RMSE en 
( 3.2.4 SAARA EUER E AM 苗 期 0.77 0.70 0.77 4.59 
T LAS A" RR ER AAEE, Hag HE 封 行 期 0.81 0.79 0.82 6.39 
© SRAZI, EMA, ZMR M oss ama 
线性 拟 合 ， 结 果 见 表 3。 株 高 监测 精度 RR 为 ”_ 人 和 网 一 — % —— 56 
0.70~0.86, RMSE 控制 在 4.59~6.86 cm iG Hi], 3 = 0.0093x + 0.3405x, + 10.312 (20) 
种 回归 模型 的 拟 合 效果 存在 差异 ， 总 体 精度 表现 = 0.0034x, + 0.5128x, + 19.105 (21) 
为 多 项 式 回 归 > 线 性 回归 > 指数 回归 。 采 用 多 项 式 y, = 0.0005x, + 0.7262x, + 21.147 (22) 


回归 模型 对 世 麻 种 质 资源 全 生育 期 内 的 株 高 数据 
进行 拟 合 ， 基 于 全 生育 期 的 芝麻 株 高 提取 精度 很 
高 ， 实 测 株 高 与 Hw 之 间 的 尼 为 0.93，RMSE 为 
5.65 cm， 表 明 利 用 无 人 机 遥感 监测 芝麻 株 高 具 
有 可 行 性 。 

将 各 生育 期 数据 按 7:3 的 比例 随机 划分 为 训 
练 集 和 验证 集 ， 并 采用 训练 集 数据 构建 不 同 生 育 
期 芝麻 株 高 校正 模型 。 其 中 ， 苗 期 、 封 行 期 、 旺 
长 期 芝麻 株 高 校正 模型 分 别 为 公式 (20) ~ (22). 


表 4 为 节 麻 株 高 校正 模型 的 校正 效果 ， 由 
表 4 可 知 ， 校 正 后 验证 集 的 精度 明显 提升 ， 校 正 
后 的 Huss 与 实测 株 高 有 较 高 拟 合 度 。 
3.2.2 ”区 麻 种 质 资源 生理 参数 监测 

LR, RF, SVM, PLSR 等 4 种 机 器 学 习 算 法 
佑 测 芝 麻生 理 参 数 结果 如 表 5 所 示 。LR 在 芝麻 
SPAD 值 估 测 中 表现 最 佳 ， 测 试 集 尼 最 高 达到 
0.66, RMSE 为 2.03; 对 含水 量 估 测 表 型 也 为 最 
f£, R?=0.79, RMSE-221, XJ FE II ZE BE Ta 
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XA RRS RIE ER 
Table 4 Effect of ramie height correction model 
校正 前 校正 后 
生育 期 
R RMSE/cm R RMSE/cm 

苗 期 0.85 5.10 0.98 1.14 
封 行 期 0.77 6.06 1.00 4.85 
旺 长 期 0.90 7.01 0.10 3.18 


精度 (R= 0.78，RMSE=2.29) ， 测 试 集 的 拟 合 
旨 标 得 到 了 提高 ， 说 明 该 模型 具有 较 强 的 泛 化 能 
力 和 稳定 性 。SVM 在 芝麻 LAI 估 测 中 表现 最 佳 ， 
训练 集 和 测试 集 的 尼 分 别 为 0.74 和 0.74，RMSE 
分 别 为 0.60 和 0.63， 训 练 集 与 测试 集 的 相似 精度 
表明 该 模型 是 可 靠 稳 定 的 。 因 此 ， 本 研究 基于 遥 
感 影 像 提取 的 芝麻 生理 参数 具有 较 高 的 准确 性 ， 
后 续 研 究 将 基于 最 佳 模型 估 测 值 对 芝麻 种 质 资源 
进行 分 析 。 
表 5 世 麻 种 质 资 源 生理 参数 监测 结果 

Table 5 Monitoring results of physiological parameters of 


ramie germplasm resources 


训练 集 测试 集 


指标 。 模型 - - 
R RMSE R? RMSE 


LR 0.69 2.05 0.66 2.09 
RF 0.93 0.97 0.59 2.30 

SPAD {Ëi 
SVM 0.61 2.28 0.66 2.09 
PLSR 0.62 2.27 0.64 2.14 
LR 0.78 2.29 0.79 2.21 
RF 0.96 0.93 0.77 2.33 

含水 量 
SVM 0.76 2.40 0.78 2.30 
PLSR — 0.77 2.36 0.79 2.20 
LR 0.77 0.56 0.70 0.67 
RF 0.96 0.24 0.71 0.66 

LAI 

SVM 0.74 0.60 0.74 0.63 
PLSR — 0.74 0.60 0.73 0.62 


3.9 PS RAH EX URBIS ES 


对 蔚 麻 种 质 资源 遥感 表 型 进行 变异 性 分 析 
(图 4)。 不 同时 期 芝麻 种 质 资源 LAI 的 变异 系数 
范围 为 18.94%~23.42%,， 平均 变 异 系 数 为 
20.83%; How 的 变异 系数 范围 为 20.06% 一 


27.27%， 平 均 变 异 系 数 为 24.61%; 株数 的 变异 
系数 为 35.48%。 以 上 结果 表明 ， 苍 麻 种 质 资源 
在 遥感 提取 的 LAI、 株 高 、 株 数 结果 上 具有 和 较 大 
差异 。 在 实际 生产 中 ， 株 高 、 有 效 株数 是 构成 芝 
RPO AY EEA), LAT BE IRA RAP AE 
结构 以 及 光合 作用 能 力 ， 也 与 最 终 产 量 密切 相 
关 。 因 此 ， 考 察 营 麻 资 源 的 株 高 、 株 数 、LAI 差 
异 将 有 助 于 蔚 麻 产量 的 评 佑 和 高 产品 种 的 选 育 。 
含水 量 的 变异 系数 范围 为 2.3% 一 4.59%， 不 同 种 
质 资 源 间 的 差异 性 随 着 生育 期 的 后 移 逐 渐 增 大 , 
在 旺 长 期 (4 月 20 日 ) 达到 最 大 值 4.95%， 这 是 
由 于 该 时 期 气温 上 升 ， 降 水 减少 ， 土 壤 含 水 量 降 
低 ， 相 比 蕊 麻生 长 早期 ,该 时 期 可 以 看 作 是 干旱 
胁迫 。 不 同 芝 麻 种 质 资源 含水 量变 异 范围 的 增 大 
表明 品种 耐 旱 能 力 存在 差异 。 


40 


W3.158 3238 3.290 4.06 4.128 4.2 


变异 系数 /% 


LAI SPAD ” 株 高 /cm KE% ”株数 / 株 
遥感 表 型 特征 
图 4 基于 天 感 衍 生性 状 的 芝麻 种 质 资源 差异 性 分 析 
Fig.4 Difference analysis of ramie germplasm resources 


based on remote sensing derived characters 


3.4. STER UL IER JR REI eB f e 


对 基于 遥感 的 芝麻 种 质 资源 表 型 进行 主 成 分 
分 析 (图 5)， 两 个 主要 成 分 解释 了 所 有 种 质 资 源 
性 状 超过 43.39% 的 总 变异 。 第 一 主 成 分 中 ， 株 
高 和 LAI 的 载荷 量 较 高 ， 两 者 都 是 反映 芝麻 结构 
特征 的 重要 指标 ， 因 此 主 成 分 1 (因子 1) 可 用 
于 蔡 麻 种 质 资源 结构 特征 评价 。 不 同时 期 LAI 和 
SPAD 值 在 第 二 主 成 分 中 具有 和 较 高 的 贡献 率 ， 
LAI 能 够 反映 芝麻 光 截 获 能 力 ， 叶 绿 素 含 量 也 与 
植物 光合 作用 密切 相关 ， 因 此 ， 主 成 分 2( 因 
子 2) 可 以 作为 高 光 效 芝麻 资源 的 筛选 指标 。 
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Fig. 5 Principal component classification of remote sensing 


phenotypes of ramie germplasm resources 


壬 芝麻 种 质 资源 多 样 性 分 析 基 础 上 ， 进 一 步 
tt sii i 
综合 性 指标 ， 根 据 综 合 性 状 能 够 获得 关于 基因 型 

异 的 补充 信息 ， ny 
坐标 系 中 的 分 布 如 图 6 所 示 。 根 据 种 质 资源 所 在 
象限 ， 可 以 判断 特定 性 状 对 应 的 最 强 影 响 。 越 靠 
近 坐 标 轴 原 点 ， 表 明 种 质 资源 在 第 一 、 二 性 状 上 
是 平衡 的 ， 而 位 于 边缘 的 种 质 资源 则 受到 特定 性 
状 的 有 影响。 如 品种 5-5 和 品种 1-9 在 因子 1 上 具有 
相似 的 性 状 , 但 是 5-5 具 有 较 低 的 因子 2 特征 值 ， 
1-9 具 有 较 高 的 因子 2 特征 值 ， 表 明 1-9 比 5-5 的 
光 效 利用 率 会 更 高 。 结 果 表 明 ， 通 过 查找 不 同 区 
麻 品 种 在 因子 1- 因子 2 坐 标 系 中 的 分 布 ， 可 以 判 
断 结构 相关 性 状 和 光 效 利用 相关 性 状 对 应 的 
影响 。 
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图 6 芝麻 种 质 资源 在 因子 1- 因 子 2 坐标 系 中 的 分 布 
Fig.6 Distribution of ramie germplasm resources in factor 1- 


factor 2 coordinate system 
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作物 种 质 资源 的 表 型 多 样 性 是 品种 选 育 的 动 
力 ， 本 研究 结合 无 人 机 遥感 技术 、 机 器 学 习 技 术 
提出 了 一 种 快速 、 无 损 的 芝麻 种 质 资源 表 型 多 样 
性 评价 方法 。 主 要 结论 如 下 : 

(1) 基于 无 人 机 遥感 的 芝麻 种 质 资源 表 型 估 
测 是 可 行 的 。 田 间 实 测 表 型 数据 与 遥感 表 型 之 间 
的 拟 合 效果 较 好 ， 基 中 芝麻 株 高 估 测 精度 为 
0.93, RMSE 为 5.65 cm; SPAD {if R 为 0.66， 
RMSE 为 2.09、 含 水 量 尼 为 0.79，RMSE 为 
2.21%, LAI AYR’ 0.74, RMSE 0.63. WHI 
人 机 遥感 获取 的 低空 高 分 辩 率 影像 信息 不 仅 能 实 
时 、 高 频 地 为 分 析 作 物 长 势 状况 提供 数据 支撑 ， 
还 能 作为 现 有 地 面 、 航 空 、 卫 星 遥 感 的 补充 ， 构 
建 更 为 全 面 、 立 体 的 农业 遥感 监测 系统 2 。 

(2) ELR, SVM, RF, PLSR 分 别 构建 
了 基于 无 人 机 和 遥感 的 芝麻 生理 参数 估 测 模型 ，4 
种 机 吉 学 习 算 法 的 估 测 效果 相当 ，LR 在 芝麻 
SPAD 值 、 仿 水量 佑 测 中 表现 最 佳 ， 而 SVM 在 
LAI 佑 测 中 表现 最 佳 。 

(3) 提取 的 遥感 表 型 具有 反映 芝麻 种 质 资源 
基因 型 差异 的 潜力 ， 其 中 LAI、Husw 和 株数 的 变 
异 系数 分 别 达 到 20.8396, 24.61% 和 35.48%， 说 
明基 于 无 人 机 遥感 技术 的 作物 种 质 资源 表 型 多 样 
性 调查 是 可 行 的 。 

(4) 对 遥感 表 型 进行 主 成 分 分 析 ， 因 子 1 可 
用 于 芝麻 种 质 资源 结 构 特 征 评价 ， 因 子 2 可 以 作 
为 高 光 效 芝麻 资源 的 筛选 指标 ， 主 成 分 分 析 结果 
为 芝麻 种 质 资源 规模 化 高 效 鉴定 ， 基 因 型 分 析 与 
育种 利用 创造 了 可 能 。 

利用 无 人 机 遥感 技术 快速 、 无 损 、 高 效 提取 
大 规模 田间 表 型 数据 对 作物 品种 筛选 以 及 加 速 作 
物 育 种 进程 具有 重要 意义 。 但 如 何 挖掘 遥感 数 
据 ， 从 中 茜 取 可 靠 的 、 有 意义 的 生物 学 知识 仍然 
是 未 来 研究 的 重点 。 在 今后 的 研究 中 ， 应 结合 特 
定 环境 下 作物 的 实际 育种 目标 ， 通 过 引入 最 新 的 
机 器 学 习 、 深 度 学 习 等 人 工 智 能 方法 提高 遥感 表 
型 的 估 测 精度 。 另 外 ， 作 物品 种 选 育 高 度 依赖 基 
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因 型 和 环境 因素 ， 今 后 研究 应 更 加 重视 表 型 分 析 
技术 与 基因 型 技术 的 融合 ， 以 期 筛选 特定 环境 胁 
迫 下 适应 性 更 强 的 作物 品种 。 
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High Quality Ramie Resource Screening Based on UAV 
Remote Sensing Phenotype Monitoring 


FU Hongyu, WANG Wei, LIAO Ao, YUE Yunkai, XU Mingzhi, WANG Ziwei, 
CHEN Jianfu, SHE Wei, CUI Guoxian 


(College of Agronomy, Hunan Agricultural University, Changsha 410128, China) 


Abstract: Ramie is an important fiber crop. Due to the shortage of land resources and the promotion of excellent varieties, the 
genetic variation and diversity of ramie decreased, which increased the need for investigation and protection of the ramie germ- 
plasm resources diversity. The crop phenotype measurement method based on UAV remote sensing can conduct frequent, rapid, 
non-destructive and accurate monitoring of different genotypes, which can fulfill the investigation of crop germplasm resources 
and screen specific and high-quality varieties. In order to realize efficient comprehensive evaluation of ramie germplasm pheno- 
type and assist in screening of dominant ramie varieties, a method for monitoring and screening ramie germplasm phenotype 
was proposed based on UAV remote sensing images. Firstly, based on UAV remote sensing images, the digital surface model 
(DSM) and orthophoto of the test area were generated by Pix4dmapper. Then, the key phenotypic parameters (plant height, 
plant number, leaf area index, leaf chlorophyll content and water content) of ramie germplasm resources were estimated. The 
subtraction method was used to extract ramie plant height based on DSM, while the target detection algorithm was applied to ex- 
tract ramie plant number based on orthographic images, and four machine learning methods were used to estimate the leaf area 
index (LATI), leaf chlorophyll content (SPAD value) and water content. Finally, according to the extracted remote sensing pheno- 
typic parameters, the genetic diversity of ramie germplasm was analyzed by using variability analysis and principal component 
analysis. The results showed that: (1) The ramie phenotype estimation based on UAV remote sensing was effective, with the fit- 
ting accuracy of plant height 0.93, and the root mean square error (RMSE) 5.654 cm. The fitting indexes of SPAD value, water 
content and LAI were 0.66, 0.79 and 0.74, respectively, and RMSE were 2.03, 2.21 and 0.63, respectively; (2) The remote sens- 
ing phenotypes of ramie germplasm were significantly different, as the coefficients of variation of LAI, plant height and plant 
number reached 20.82%, 24.61% and 35.48%, respectively; (3) Principal component analysis was used to cluster the remote 
sensing phenotypes into factor 1 (plant height and LAI) and factor 2 (LAI and SPAD value), factor 1 can be used to evaluate 
the structural characteristics of ramie germplasm resources, and factor 2 can be used as the screening index of high-light effi- 
ciency ramie resources. This study could provide references for crop germplasm phenotypic monitoring and breeding correla- 
tion analysis. 


Key words: ramie; diversity of germplasm resources; phenotype; UAV remote sensing; digital surface model; machine learning 
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